Convection in a Very Compressible Fluid: Comparison of Simulations With 

Experiments 

H. Meyer 1 , A. Furukawa 2 , A. Onuki 2 and A.B. Kogan 1 ’* 

1. Department of Physics, Duke University, Durham, NC 27708-0305, USA 
2. Department of Physics, Kyoto University, Kyoto 606-8502, Japan . 

*Present address: Department of Physics, MIT, Cambridge, MA 02139-4307, USA. 


The time profile A T(t) of the temperature difference, measured across a very compressible fluid 
layer of supercritical 3 He after the start of a heat flow, shows a damped oscillatory behavior before 
steady state convection is reached. The results for A T(t) obtained from numerical simulations 
and from laboratory experiments are compared over a temperature range where the compressibility 
varies by a factor of « 40. First the steady-state convective heat current j conv as a function of the 
Rayleigh number Ra is presented, and the agreement is found to be good. Second, the shape of 
the time profile and two characteristic times in the transient part of AT(f) from simulations and 
experiments are compared, namely 1) t 0 sc, the oscillatory period and 2) t p , the time of the first peak 
after starting the heat flow. These times, scaled by the diffusive time td versus Ra, are presented. 
The agreement is good for tosc/ri), where the results collapse on a single curve showing a power- 
law behavior. The simulation hence confirms the universal scaling behavior found experimentally. 
However for tp/m, where the experimental data also collapse on a single curve, the simulation results 
show systematic departures from such a behavior. A possible reason for some of the disagreements, 
both in the time profile and in t p is discussed. In the Appendix a third characteristic time, f m , 
between the first peak and the first oscillation minimum is plotted and a comparison between the 
results of experiments and simulations is made. 


I. INTRODUCTION 

In a Raleigh-Benard (RB) cell, the start of a constant heat flow across the fluid layer produces an initial rise in the 
temperature difference across this layer, A T(t), with a transient profile determined by the fluid convection dynamics, 
and which then tends to a steady state value, labeled AT. Recently such heat flow experiments were carried out [1,2] 
on a very compressible fluid in its convective state, supercritical 3 He, along the critical isochore < p >= p c , where 
the critical temperature is T c = 3.318 K. The fluid layer height in the experiments was T=0.106 cm and the diameter 
of the cylindrical cell was 5.7cm. For this large aspect ratio, the predicted critical Rayleigh number is i?a c =1708, as 
was confirmed by the experiments [1]. Over the reduced temperature range 0.009< e = (T — T c )/T c < 0.2 where the 
experiments were analyzed, the isothermal compressibility increases by a factor of « 40 as T c is approached. The 
substantial change in the fluid properties along p c is reflected in large changes of the transient profile, where damped 
oscillations were observed after the first peak of A T(t) for e > 0.009. 

This paper compares the results of experimental data, and of simulations extending those of refs. [3,4], for the time 
profile A T(t) in the regime where the fluid approaches steady-state convection. Reference is also made of recently 
published simulation results by Amiroudine and Zappoli [5]. In the simulations done in refs. [3,4], two new terms 
are added in the heat conduction equation; the first takes into account the adiabatic heating taking place throughout 
the cell (the “piston effect”), and the second accounts for the adiabatic temperature gradient effect within plumes 
leading to the Schwarzschild criterion of the convection onset in compressible fluids [6]. In the work of ref. [5] the 
Navier-Stokes (NS) equation, together with an enthalpy equation were used, without a specific term to simulate the 
piston effect, since the latter is implicitly taken into account in the work of the pressure forces term in the enthalpy 
equation. The simulations were made in two dimensions, where the fluid was contained in a cell with the same height 
L as in the experiment. The aspect ratio was 4 in refs. [3,4] and 2 in ref. [5] with the periodic boundary condition in 
the horizontal direction in order to reproduce the conditions of the experimental cell with a large aspect ratio. The 
simulations for e = 0.05 in ref [4] extended until [Ra — Ra c ] k 4 x 10 6 and until ss 1.7 x 10 5 for other values of e. 
Simulation results in this paper, except those of e = 0.05 in refs. [3,4], are newly obtained using the scheme in ref. [4]. 

The outline for the remainder of this paper is as follows: First, a general discussion of the profiles A T(t) will be 
made, with presentation of some examples and an appraisal of the degree of agreement between experiments and 
simulations. Second, a comparison of the steady state results from experiments and simulations will be presented, 
expressed in terms of the convection heat current j conv versus the reduced Rayleigh number ra* = [Ra — Ra c ]/Ra c . 
Third, a comparison of the transients from experiments and simulations will be made by the examination of two 
“characteristic times”, the time t p of the first peak after the start of the heat current, and t osc , the oscillation 
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period in the damped oscillatory decay. Both will be presented in a scaled form, again versus [Ra — Ra c \. This 
quantitative comparison illustrates well the agreements and deviations between simulations and experiments, which 
will be discussed in the Conclusions. In the Appendix, a third characteristic time, t m , between the first peak and the 
first minimum in the A T(t) profile is shown for both the experiments and the simulations. 


II. GENERAL OBSERVATIONS ON THE PROFILES AT VERSUS TIME. 

We briefly emphasize, as was done in refs [1,2], that the experiments were done under conditions where the strat- 
ification from gravity was small and where the temperature changes across the fluid were kept small enough that 
the changes in the fluid properties across the fluid layer were smaller than a few percent. Thus the conditions for 
the approximations in a Boussinesq-like fluid in the momentum equation were maintained. In Fig. 9 in ref. [1] a 
representative evolution is shown of the observed profile A T(t) at e = 0.05 as a function of the heat current q. Both 
experiments and simulations show over a certain range of heat flow Q and of temperature a damped oscillatory profile 
A T[t) as shown in Figs. 1 and 2 of ref. [3], Figs. 1 of ref. [4], both at e = 0.05 and Fig. 2 of ref. [5] at e = 0.01, which 
are not reproduced here. Considering that the experimental data have not been corrected for the time lag introduced 
by the temperature recording instrumentation with a time constant of r = 1.3 s., the qualitative agreement is quite 
satisfactory. However as e increases and the compressibility decreases, the disagreement between the transient regime 
of experiments and simulations becomes sizeable. This can be seen in Fig.l at e = 0.2 for two values of the heat flow 
q « 3.7xl0~ 7 and 2.1xl0 -7 W/cm 2 . By contrast the steady-state value for 



t(s) 

FIG. 1. Plots of the profile A T(t) versus time for two values of the heat flow q (in 10 -7 W/cm 2 ) at e = 0.2 and comparison 
between experiments and simulations a) q = 3.89 (expt), 3.67 (simul.) and b) q = 2.16 (expt), 2.10 (simul.). 

AT reached in both the experiments and in simulations for the same value of q remains closely the same. In the 
experimental trace for the lower value of q , no damped oscillations are seen, but rather a non-exponential decay of 
the overshoot. This is the regime labeled “truncated oscillations”, described in Fig. 2 of ref. [2]. (See also Appendix 
A). Simulations, however, always show damped oscillations in the convective state. 
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FIG. 2. A representative profile A T(t) after starting the heat current, with the definition of the times f p , t m and t Q sc and 
the slow exponential relaxation to the steady-state value of AT (dashed line) with a time constant Ttaii • The recording is for e 
= 0.02 with the heat flow q = 1.69 10 ' W/cm 2 ). 

In Fig. 2, we present an enlarged segment near the peak of AT(t') for a representative trace at e = 0.02, which exhibits 
and defines various characteristic times in the transient. We will consider two of these times in the main part of this 
paper, namely 1) the peak t p of the first peak after the start of the heat current and 2) the period f osc of the damped 
oscillations. 

Furthermore there is t rn , the interval between the first peak and the first minimum in the regime of damped 
oscillations or between the peak and the “kink” to the flat portion in the regime of “truncated oscillations”. The 
discussion of t m is deferred to Appendix B, because it is more complex than for t osc and t p . Finally, the profile tends 
asymptotically from below to the steady state value AT with a relaxation time Ttaii, obtained by fitting this transient 
tail with a simple exponential, and we note that this feature has not been detected in the simulations either in ref. 
[4] or in ref. [5]. 

The transient shape, principally at short times where the changes in A T(t) are rapid, is affected by the instrumen- 
tation time constant, and also the modulated minimum shown in Fig. 2 makes an analysis of the oscillation amplitude 
decay rate uncertain. However inspection of the many recorded experimental traces shows that for a given e (or 
compressibility), the rates for both the oscillations and their amplitude decay increase with the heat current. For 
values large enough of q at a given e, the first peak and the oscillations become attenuated and averaged out. (See Fig. 
4 of ref [2]). The likely cause for this observation is that t osc becomes comparable or smaller than the instrumentation 
time constant. An averaging effect of the oscillations due to a negative interference between non-synchronous plumes 
that results from the large lateral dimension of the experimental cell was suggested in ref. [4] . This suggestion might 
be very relevant in the regime e < 0.009 where no oscillations could be observed, even when their expected period 
was well above that of the instrumental time constant (See ref. [2], Section III). 


III. COMPARISON OF STEADY-STATE RESULTS IN EXPERIMENTS AND SIMULATIONS 

In ref. [1] the steady-state experimental results were presented in terms of the dimensionless convective heat current 
jconv versus ra *_ n ere j conv j g ^e ra fj 0 0 f convective portion of the heat current to that conducted through the 
fluid at the transition to convection, which leads to the relation [7] 

j conv = (JVu-l)(ro* + l) (1) 

where Nu is the Nusselt number. Along the critical isochore of a fluid, where the compressibility diverges as T c 
is approached, Gitterman and Steinberg [8] have shown that for the fluid onset of mechanical instability there is a 
crossover from the Rayleigh - to the Schwarzschild (or “ adiabatic temperature gradient”) criterion as T c is approached. 
Then it has been shown ( [9] and references therein) that at the onset of convection, the temperature drop across the 
fluid layer is given by 
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AT onset — A Tp + A T a d (2) 

where ATr = Ra c x vDp /apgL 3 , the familiar Rayleigh term, and A T a d = LgTap/Cp, the adiabatic temperature 
difference contribution (for ATr and A T a d see for instance Tritton’s monograph [10]). Here v is the kinematic 
viscosity, Dp the thermal diffusivity, ap the isobaric thermal expansion coefficient, g the gravity accelation, L the 
height of the fluid layer and Cp the specific heat at constant pressure. 

As has been described in some detail in refs. [1,2], Nu and Ra have to be corrected for the contribution from the 
adiabatic temperature gradient. One then obtains 


Ra 


corr 


Ra(AT - A T ad ) 
AT 


and Nu corr 


( AT di ff — A T ad ) 
(AT - A T ad ) 


( 3 ) 


and hence 


ra* COII = ( Ra c 


r/Ra c - 1) and = ( Nu COII - l)(ra* orr + 1) 


( 4 ) 


Here ATdiff is the temperature drop across the fluid in the diffusive regime for the same heat current producing the 
observed AT. Because both and ra* orl . vary over about five orders of magnitude for the range covered by the 
data, a more sensitive way is to present the ratio jco™ /ra* orl versus ra* orr . Furukawa and Onuki [4] theoretically 
justified the validity of the scaling relations in terms of these corrected quantities. This is done in Fig. 3a where the 
points recorded close to the transition to convection and showing rounding instead of a sharp convection onset, have 
been omitted. Furthermore the data for e <0.009, where no damped oscillations were obtained in the transients, and 
for which no simulations were carried out, have not been used in this figure. As can be seen, within the scatter all 
the data points nearly collapse on a common curve. For ra* orr <1, the data extrapolate to a horizontal line with an 
amplitude of 1.3 ± 0.1. This asymptotic result, = 1.3ra* orr , which represents data slightly above the onset of 
convection, has been discussed in ref. [1] where it had been concluded that the amplitude is consistent with a straight 
roll convection prediction [11]. In Fig. 3b, the experimental data are replaced by a solid curve representing its average, 
and the results from simulations are shown by the symbols at various values of e. In the intermediate range of ra* orr 
these data points collapse on a curve slightly above the experimental average, and also within their scatter they tend 
to the same limiting amplitude of 1.3, as do the experiments. Therefore they also imply straight roll convection, a 
picture which is plausible in the 2D simulation which represents a cross-section of rolls in a geometry with parallel 
vertical periodic boundaries. 
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FIG. 3. The convective heat current j divided by ra* orr =[Ra COII — Ra c ]/Ra c , versus ra* orr , both corrected for the 
adiabatic temperature gradient contribution, a) Experimental data at various values of e. b) Data from simulations at different 
values of e shown by symbols. The solid line is the average of the experimental data in a). 

IV. TRANSIENT CHARACTERISTIC TIMES IN EXPERIMENTS AND SIMULATIONS 


A. Background information 


The mechanism responsible for the damped oscillations for a highly compressible fluid at constant average density 
has been discussed most recently by Furukawa and Onuki [4] and by Amiroudine and Zappoli [5]. From their 
simulations, these authors presented a detailed analysis of the time evolution of the temperature in various locations 
in the fluid layer. We refer to their description on how maxima and minima of A T(t) are produced by the “piston 
effect” that leads to the vertical flow of successive “warm” and “cold” fluid masses. The snapshots of the simulations 
by Chiwata and Onuki [3] had already pictured the formation of plumes during these processes, and the evidence of 
warm and cold vertical flows. 

In the analysis of the experimental data [2] , two characteristic times describing the remarkable oscillatory behavior, 
t rn and t osc , were discussed, which are shown in Fig. 2 for a representative A T(t) profile. The experimental results for 
the time t v have not been analyzed before and are presented here for the first time. In order to keep the analysis of 
the results tractable, their discussion will not include t m , which will be dealt with in Appendix A. In the experimental 
data analysis, the relaxation time T ta u to the steady state convection was also presented in ref. [2] . All the data of t osc 
and T tai i at the various reduced temperatures e, scaled by the diffusion time to = L 2 /AD T , and plotted versus the 
Rayleigh number difference [Ra corr — i?a c ] were found to collapse within a scatter of ±15% on two respective curves, 
an unexpected result. 


e 

C P /C V = 7 

III 

to 

t d (s.) 

0.01 

119 

1.01 

265 

0.02 

57.7 

1.03 

134 

0.03 

38.1 

1.06 

92.0 

0.04 

28.5 

1.08 

69.9 

0.05 

22.8 

1.10 

57.0 

0.07 

16.5 

1.13 

42.2 

0.10 

11.8 

1.19 

31.8 

0.20 

6.48 

1.35 

19.8 


Table 1. The ratio C p /C v = 7 for 3 He along its critical isochore for several values of the reduced temperature e, 
the ratio B(j) of the diffusive times, calculated from Fig.2 ref.12 (Behringer et al.) and the diffusive time td- 

We now discuss the choice of the scaling time in the data for t osc and t p versus the Rayleigh number. As mentioned 
above, the diffusion time tn was used as the scale for t osc in the plot versus [Ra COII — Ra c ) . This time to was first used 
in the Navier Stokes equation leading to the simulations of ref. [3] under the conditions of constant average density 
and valid for C p /C v = 7 >> 1. In the absence of convection, to can be related to the relaxation time tq = L 2 /tt 2 Dt 
of the lowest diffusion mode in a RB - or in a standard thermal conductivity cell [12]. However over the experimental 
and simulation range where 0.2 > e > 0.009, the condition 7 ^ 1 is only progressively realised as e decreases to 
0.01. In ref. [12] expressions for the solution of To as a function of 7 were derived. As 7 — > 1, To — > AL 2 /tt 2 Dt, 
and this last value is the same as for a fluid relaxing at constant pressure. In general to ( 7) = B(j)L 2 /V 2 Dt, where 
Bin) = ( 7r / ( ?o-b) 2 , obtained from Fig.2 of ref [12] with go the wave number of the lowest mode. B{ 7) is presented 
in Table 1 for several values of e relevant to the experiments and simulations in this paper. As can be seen in this 
Table, the ratio to( 7 )/t 0 (7 = 00) decreases with e and tends to 1 for e <0.01. The use of the time td = to( 7)7 t 2 /4 as 
a scale for the convection transient characteristic times is therefore preferable to tu- It turns out that the choice of 
td as the scaling time improves the collapsing of the experimental data points for both t osc and t p . In the following 
subsections we present and compare these times obtained experimentally and from simulations. 


B. the oscillatory period t 0 sc 

In ref. [2], Fig. 4 presented t osc scaled by t D and plotted versus [Ra COII — Ra c \. The substitution of t d as a scale 
gives a better collapsing of the data points than does tp, particularly at the lower values of [Ra COI r — Ra c ] where most 
of the data points at the higher values of e lie. Within the scatter of ±15% the data for 10 3 < [f?a CO rr — Ra c ] < 5 x 10 5 
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can be represented by a power law with an exponent of -0.52 ± 0.02. In Fig. 4a we present these data and in Fig. 4b 
its average, obtained by a fit to a power law with a corrective term for the larger Ra values. Also in Fig. 4b we 
show by symbols the results from the simulations by the Kyoto group (two of the present authors, A.F. and A.O. and 
Chiwata) and by Amiroudine and Zappoli [5] at various values of e. The agreement appears to be good, though the 
simulations indicate some small systematic deviations from collapse on a single “universal” line. The basic reason for 
the apparent power law with an exponent of -1/2 remains to be understood. 




(Ra° u - Ra c ) 

FIG. 4. The oscillation period t osc and the time Ttaii for the relaxation to the steady-state, both scaled by the diffusion time 
td, versus [Ra corr — Aa c ]. a): Experimental data at various values of e. b): Symbols show the simulation data at different 
values of e, and the solid line is the average of the experimental data in a). 


C. The location of the first peak at f p 

In the experiments, the initial rise of the measured transient A T(t) after the start of the heat flow across the 
fluid layer is affected by the time constant of the thermometer circuitry, r = 1.3 s. as described in ref. [1]. This is 
especially so at short times, when the A T[t) increases rapidly, and as a result the measured value of the time at the 
first peak has to be corrected. This was done by comparing the shift in time between the calculated rise of A T(t) in 
the conducting regime (Eq.3.3 of ref. [4]) and the recorded curve for several values of e and q. This shift was between 
2 and 3 s. and t p was obtained after a crude correction was made by subtracting 5t = 2 s. from the measured time 
at the peak. The times t p used in this analysis ranged from « 100 to « 7 s. Fig. 5a shows the scaled representation 
tp/T~D of the experimental times versus [f?a CO rr — Ra c ], There is excellent collapse of the data for all the values of e, 
even extending to the region of e < 0.009 where no damped oscillations are observed [2]. In Fig. 5b the solid curve 
represents the average of the experimental data, as obtained by a polynomial fit. Symbols show the results from 
simulations at various values of e. The data from ref. [3] at e = 0.05 
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FIG. 5. The time of the first peak i p , scaled by the diffusion time td, versus [f?a CO rr — i?a c ]. a): The experimental data, 
where the symbols denote the various values of e at which the experiments were carried out. b): The curve representing the 
average of the experimental data from 5 a), compared with the results from simulations, shown as symbols at the various values 
of e. 


and those in the present research are found to be internally consistent. However, contrary to the experimental results, 
the simulation data do not collapse on a single curve. This discrepancy is expected from the profiles A T{t) shown 
in fig. 1 where the peak in the simulations lies at substantially longer times t p than for the experiments. A possible 
source for the discrepancy between simulations and experiments will be considered below. 


V. DISCUSSION 

First we present general comments on transient observations in Rayleigh-Benarcl convection. After the start of a 
heat current at constant pressure, an overshoot in A T(t) is routinely observed in Boussinesq fluids. The origin of the 
overshoot is a certain “inertia” of the fluid immediately after the heat flux is applied: Even though the final state 
is convection, initially the fluid remains stationary and A T(t) follows the solution of the thermal diffusion equation. 
As the fluid begins to move, the trace shows an overshoot with a peak at t v followed by a transition to the steady 
state value in the convective regime as described, for example, by Behringer [13]. Thus, t p approximately describes 
the time that it takes for the layer to develop convective motion. 

In a compressible fluid, the non-convecting state during the initial transient is expected to persist up to a value 
AXlnstab greatly exceeding the stability threshold AT ons given by Eq.2, derived for a linear temperature and density 
distribution in the vertical direction. This result for AT| ns tab was obtained by El Khouri and Carles [14] via a 
linear stability analysis and is a direct consequence of the strongly inhomogeneous vertical density and temperature 
distribution with pronounced boundary layers generated by the Piston effect [15]. A clear example can be seen in 
Fig. 1(a) of ref. [4], which presents transient curves for e = 0.05. The simulated trace essentially coincides with the 
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zero-gravity result (Eq.3.3 in ref [4]) up to A T(t) of order 280 [iK, while the stability criterion (Eq.2 in the present 
paper) gives only 7 /iK. Another example is obtained by comparing the simulations in Fig.2 of ref [5] at e = 0.01 with 
calculations under zero gravity, which also shows the two curves coinciding until close to t p . Hence the simulations 
[4,5] are qualitatively consistent with the predictions of linear stability analysis. 

We now recapitulate the main findings of the comparison between experimental data and simulations for a su- 
percritical very compressible fluid, 3 He. Starting with t osc , which is the characteristic time determined farthest away 
from the start of the heat flow, both experiments and simulations in the scaled representation are in good agreement. 
By contrast, as shown in Fig. 5, t p is systematically larger in the simulations than in the experiments, the difference 
becoming more important as e increases. An intriguing puzzle is why all the experimental t p data at the various 
values of e can be cast into a scaling representation, while the simulation data cannot. 

This disagreement for t p is surprising: at large e where the compressibility of the fluid has become smaller, a 
“simple” Boussinesq behavior should be recovered. A possible origin of this discrepancy is the imperfection of the 
temperature control of the top plate during the experiment. The signal from the control thermometer [1] is recorded 
by a circuit with a time constant of 1.25 seconds and 6 dB/octave rolloff. Therefore, temperature fluctuations with 
frequencies above a few Hz cannot be detected. The average temperature noise is estimated to be about 1 [iK/\jHz 
r.m.s. Because of very high thermal conductivity of the top plate material (OFHC copper) it seems reasonable 
to assume that this perturbation produces no horizontal temperature gradients. It is therefore possible that small 
parasitic fluctuations of the top plate temperature could speed up the development of the convecting state, hence 
producing an overshoot with a smaller t v than the simulations do, which have no noise. We would anticipate that once 
the convection is almost fully developed, the influence of the fluctuations mentioned above would not be noticeable 
any more. Hence they might not affect the period t osc of the damped oscillations, which would explain the good 
agreement between experiment and simulations and also that in the steady-state condition [16]. 

We mention here that in the experiments [1] a rounding of the onset point on the steady-state A T{q) measurements 
was reported and also attributed to the top plate noise. The rounding was found to become more and more pronounced 
as the critical point point was approached, i.e. the effect of the fluctuations would have to be increasing with decreasing 
e. The effect of time-dependent boundary conditions on convection in a compressible fluid appears to be an interesting 
and open question. We suggest future studies, both experimental and numerical ones, that would focus on the transient 
response of a fluid layern a R-B cell to an externally imposed perturbation of the top surface temperature. 


VI. CONCLUSIONS 

A systematic comparison of the results from experimental and 2D simulation convection studies of a compressible 
fluid has been reported. This fluid is supercritical He 3 along the critical isochore and over a temperature range 
where the compressibility varies by a factor of ss 40. The temperature profile A T(t) across the fluid layer in a 
Rayleigh-Benard cell after the start of a heat flow was investigated in both experiments and simulations. The damped 
oscillations in the transient after the start of the heat flow, and with a period t osc are of particuliar interest. They 
are a consequence of vertical mass flows that result from the “piston effect” triggered by the plumes - both moving 
up and down. 

The comparison of the steady-state results, expressed in terms of the convection current versus the Rayleigh 
number, shows good agreement in general. However, the comparison of the transient results shows some systematic 
discrepancies, which appear to become more important as the compressibility decreases (i.e as the distance from the 
critical point increases). This can be seen clearly by examining the respective A T(t) profiles at various temperatures. 
The agreement is best for the oscillation period where the simulation results and the experimental data can be 
represented in scaled form versus the Rayleigh number. The absence of noise in the simulations, in contrast to a 
physical system, might be a possibility for the discrepancy. In addition, as mentioned earlier in this paper, the 
simulations carried out independently by two research groups [4,5] do not detect the slow relaxation of A T(t) to the 
steady state, or the region of “truncated oscillations” observed in the experiments. These discrepancies remain to be 
understood. 
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VIII. APPENDIX 


A. Boundary between regimes with damped and “truncated” oscillations. 

It is useful to present this boundary in terms of the Rayleigh number difference [Raj?° r 'J nd — Ra c ] versus e, where 
the transformation from the former representation A T\, oun( i(e) to the present one was done via Eq.6 of ref. [2], The 
zone diagram in this format shows that the boundary height remains constant at [Ra h c °^ ld — Ra c ] = 700 ±200 for 
0.05 < e < 0.2, and then rises steeply as e further decreases. This is shown in Fig. 7. 



10' 2 10' 1 


oo 

FIG. 6. The boundary between the regimes showing damped oscillations (top region) and truncated oscillations (bottom 
region), expressed in terms of [Ra COII — Ra c ] versus e. 

B. The time f m in the profile A T(t) 

The discussion of t m , the time defined in Fig. 2, is more complex than for t osc and t p . In Fig. 7, both the results 
for t m / t d versus [Ra COII — Raj from the analysis of the experimental and the simulated transient profiles A T(t) 
are presented. Both sets of results are clearly separated by a roughly constant amount on the vertical logarithmic 
scale. Hence we can use the same symbols for the values of e for the experiments and the simulation. For the sake 
of simplicity, we do not distinguish between the various authors in this figure. We note, however, that the simulation 
points for [Ra cor r ~ Raj = 5.96 xlO 4 and 4.78 xlO 5 at e = 0.01 are those by Amiroudine [5]. Here, in contrast with 
tosc-, the tm data from experiments show a small but systematic departure from collapsing on a single curve at the 
lower values of [Ra corr — Ra c \ . However the scaled results from simulations at the various values of e collapse on one 
curve within the scatter. The trend is the same for both experimental and simulation data : the limiting slopes of 
the curves at the extremities of [Ra corr — Ra c ] correspond to power laws with exponents of approximately n = -1 and 
n = -0.5, respectively, at the low - and at the high end of the Ra numbers. This is shown in Fig. 7. 

We note that 1) the experimental data points for [Ra corr — Ra c } below the values of the boundary shown in Fig. 
6 are in the regime of “truncated” oscillations. They form a smooth continuation of the curves at higher values of 
[Ra corr — I?a c ], namely in the regime of damped oscillations. 2) The experimental data for e > 0.07, which collapse on 
a single curve are those where the boundary in Fig. 6 is at [Ra corr — RaJ« const = 700. The departure from collapse 
from the first onto a second curve takes place when the boundary in Fig. 6 begins to rise steeply for e < 0.07. 3) It is 
intriguing that the logarithmic vertical separation of the experimental and simulation curves corresponds to roughly 
a factor of 2. This means that for the same value of [ Ra corr — Ra c ], the observed t m has roughly twice the value of 
t m from simulations. 

The analysis of the experimental t rn data for the limit of low values of [Ra corr — Ra c ] and the comparison with a 
simple predicted expression by Chiwata and Onuki [3] was already presented in Figs. 7 and 8 of ref. [2]. The exponent 
of -1 found both for the experiments and the simulations is consistent with the simple prediction. 
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Ra corr -Ra c 

FIG. 7. The time t m between the first peak and the first minimum, scaled by the diffusion time td, versus [7?a C orr — Ra c ], 
The experimental data are located on the upper two almost overlapping trajectories of symbols. Those from the simulation, 
showing a collapse, are located on the lower trajectory of symbols. The symbols denote the various values of e at which both 
experiments and simulations were carried out. 


[1] A.B. Kogan and H. Meyer, Phys. Rev. E 63, 056310 (2001). 

[2] H. Meyer and A.B. Kogan, Phys. Rev. E 66,056310 (2002). 

[3] Y. Chiwata and A. Onuki, Phys. Rev. Lett. 87, 144301 (2001). 

[4] A. Furukawa and A. Onuki Phys. Rev. E 66, 016302 (2002). 

[5] S. Amiroudine and B. Zappoli, Phys. Rev. Lett. 90, 105303 (2003), and S. Amiroudine (private communication). 

[6] L.D. Landau and E.M. Lifshitz, Course of Theoretical Physics: Vol.6 Fluid Mechanics (Pergamon, Oxford, 1959) 

[7] G. Aiders, M. Cross, P. Hohenberg and S. Safran. J. Fluid Mech.110, 297 (1982). 

[8] M. Gitterman and V. Steinberg, J. Appl. Math. Mech. USSR 34, 305 (1971). M. Gitterman, Rev. Mod. Phys. 50, 85 
(1978). 

[9] P. Carles and B. Ugurtas, Physica D 126, 69 (1999). 

[10] D.J. Tritton, Physical Fluid Dynamics (Oxford Science, Oxford, 1988), Section 14.6. 

[11] A. Schlueter, D. Lortz and F. Busse, J. Fluid Mech. 23, 129 (1965). 

[12] R.P. Behringer, A. Onuki and H. Meyer, J. Low Temp. Phys. 81, 71, (1990). 

[13] R.P. Behringer, Rev. Mod. Phys. 57, 657 (1985). 

[14] L. El Khouri and P. Carles, Phys. Rev. E 66, 066309 (2002). 

[15] These authors have also estimated ATlnstab from an analysis of experimental data corrected for the instrumental time lag, 
and have found them to be consistent with their calculations for scenarios (a) and (b) where the initial rise of A T(t) is 
not too rapid and the correction can be carried out (L. El Khouri and P. Carles, to be published, private communication). 

[16] The absence of noise in simulations, its presence in a physical system, and the resulting impact on the growth of pertur- 
bations in the convective state was pointed out to one of us (HM) by P. Carles (private communication). 


10 

25 




